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Abstract 

The standard infinite-volume definition of connected correlation 
function and particle mass in the 3-state Potts model can be im- 
plemented in Monte Carlo simulations by using C-periodic spatial 
boundary conditions. This avoids both the breaking of translation 
invariance (cold wall b.c.) and the phase-dependent and thus possibly 
biased evaluation of data (periodic b.c). The numerical feasibility of 
the standard definitions is demonstrated by sample computations on 
a 24 x 24 x 48 lattice. 
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The 3-dimensional 3-state Potts model has been extensively studied in 
particle physics context because of its relation to the finite-temperature SU(3) 
Yang-Mills theory |rj. Both have phase transitions which are characterised 
by spontaneous breaking of a global Z(3)-symmetry. A controversial issue 
for the latter has been the order of the deconfinement phase transition [Q. [| 
While two coexisting states, typical of a first order phase transition, were 
reported in Ref. ||, long-range correlations reminiscent of a second order 
phase transition were observed in Ref. |2| . The resolution of this controversy 
was first proposed in the context of the Potts model in Ref. [§j] . Using finite 
size scaling theory for bulk quantities, a first order phase transition was 
established although the correlation length displayed identical behaviour to 
that seen in Ref. |2J]. It was argued that the diverging correlation length 
was not due to a massless physical particle, but rather due to the finite-size 
effect of vacuum tunnelling. Subsequently, simulations of SU (3) gauge theory 
at finite temperature confirmed p] such an interpretation. The problem of 
determination of the physical correlation length was, however, left almost 
unsolved. Since the physical picture of a phase transition appears intuitively 
much clearer in terms of the physical correlation length, it appears desirable 
to devise a method of separating it from the tunnelling correlation length. 
Indeed, various physics applications, such as the Early Universe or heavy ion 
collsions, need a precise value of the physical correlation length in the critical 
region. 

Most simulations employ periodic spatial boundary conditions to min- 
imise surface effects. These, however, give rise to a rather complicated low- 
energy level structure in the broken phase of a theory with a spontaneously 
broken discrete symmetry. For both the Potts model and the finite tem- 
perature SU(3) theory, one has three ordered vacua and tunnelling occurs 
between them in a finite volume. Close to the transition temperature there 
is order-disorder tunnelling as well. Thus, if the mass of the physical particle 
were small, it would be hard to separate it from the multitude of tunnelling 
levels. Indeed, various prescriptions ||, |6], 0] have been used to extract the 
physical mass in the Potts model. These are, however, not satisfactory since 
the results can be sensitive to various assumptions and parameters. Thus, 
e.g., Ref. || claimed a discontinuity in the physical mass which was shown 
to decrease with increasing volume in Ref. |7|]. We therefore agree with the 
authors of Ref. |J that a useful strategy is to single out one of the three 
ordered ground states by imposing a suitable boundary condition, thus sim- 
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plifying the low-energy level structure. Our boundary condition will emulate 
spatial infinity in two respects: breaking the global Z(3) symmetry (without 
an external field) while preserving translational invariance. 

A Z(3)-symmetry breaking spatial boundary condition already exten- 
sively studied for SU(3) Yang-Mills theory is the cold wall J2|: all links lying 
in the wall are set to their trivial values. The Potts model analogue is to 
put s w = 1 for all spins on the wall. The cold wall, however, breaks the 
translational invariance. In order to reduce surface effects, one is forced to 
obtain results away from the wall where one sees little difference from the 
periodic boundary conditions. 

The virtue of periodic b.c, on the other hand, is their preservation of 
translational invariance, and avoidance of any surface interactions. With a 
cold wall b.c. the very concept of particle mass is only an approximation. 
Furthermore, particles with any kind of charge will induce anticharges on 
the wall and will thus incur surface interactions. 

Recently, C-periodic spatial boundary conditions have been developed for 
Monte Carlo simulations of (generalized) non-zero charges on finite lattices 
[§, [H], p]| . A numerical study of the SU(3) Polyakov loop using those b.c. has 



already been carried out flT|- In this Letter, we propose C-periodic spatial 
boundary conditions for extracting the physical mass of the 3-dimensional 
3-state Potts model in Monte Carlo simulations. As translational invariance 
is preserved, the concept of physical mass is clearly defined, and a precise 
determination by standard methods of statistical mechanics is possible. Con- 
nected correlation functions can be defined without any ad-hoc prescriptions 
(such as phase separation or arbitrary subtraction). 

The Hamiltonian of the Potts model with nearest-neighbour couplings is 

H = - <Wfc 

(hk) 

where cij G {0, 1, 2} defines the spin on site i. The partition function is 

M 

Link by link, H is invariant under the charge- conjugation operation 
s s where s = e 27Tia/3 and s = e ~ 2nia/3 . 
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In terms of the as, 

0^0 1^2 2-^1. 
We consider C-periodic spin configurations on an L 3 lattice, 

0"i+L = &i , (1) 

where i + L denotes the site i translated by the lattice length L in a spatial 
direction. Such configurations are conveniently visualized as a chess-board 
array of copies and C-conjugate copies of the L 3 configuration stored in the 
computer. Under translation Oi — > o~j + i of a C-periodic o configuration the 
Hamiltonian is invariant. 

The ground state of the ordered phase is characterized by <7j = on all 
sites % and is C- invariant. It is the same vacuum as with periodic boundary 
conditions and has the same energy. There are no other C-periodic vacua: if 
G{ = 1 for some i then <jj = 2 must also occur, so the minimum excitation 
energy is that of six (ij) bonds which corresponds to a particle excitation. 
A configuration with all <jj = 1 on the L 3 lattice will have all cr, = 2 on the 
adjacent C-conjugate lattices and will have the excitation energy of n^L 2 
bonds (n c bc = number of C-periodic directions). It is a many-particle state 
rather than another vacuum. 

The infinite- volume definition of the connected spin-spin correlation func- 
tion 

C(x -y) = (s x s y ) - (s x )(s y ) 

makes sense with C-periodic boundary conditions, too, since translation in- 
variance is preserved. Due to C-conjugation invariance of the Hamiltonian, 
the correlation function is real. Furthermore, 

C = R + I 
where R(x — y) = (Re s^Re s y ) — (Re s x ) (Re s y ) 
and I(x — y) = (Imsjms^) 

In any coordinate direction with a C-periodic boundary condition, Re s is 
periodic while Ims is antiperiodic |TT|] . Thus the R and I parts of the 
correlation function C are clearly distinguished here. By contrast, both R 
and / would be periodic functions if purely periodic b.c. were imposed. 
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The Fourier decomposition of C(x — y) can be interpreted in terms of 
"particle" excitations of the two-dimensional quantum field theory under- 
lying the three-dimensional Potts model. This requires that a direction of 
Euclidean time is singled out in which the boundary condition is periodic 
(finite temperature). Of the two remaining spatial directions, at least one 
should be C-periodic for lifting the degeneracy of the ordered vacua. 

For particles generated by Re s the projection onto zero spatial momen- 
tum can be implemented as usual with periodic b.c. The masses of such 
particles are best obtained from the correlation function of Res averaged 
over the two-dimensional time slices. 

For particles generated by Ims, there is a minimal spatial momentum 
of tt/L for each C-periodic direction [[J. Hence, averaging Ims over a time 
slice would mean to project onto an antiperiodic continuation of the constant 
function — an oscillatory step function with a strong contamination of higher 
momenta. The closest approximation to zero momentum in the antiperiodic 
case is a projection onto functions of minimal spatial momentum: cosx^ 
and sinx-£ for C-periodicity in the x direction, and analogously for the y 
direction. As a consequence, the mass m_ of such a particle is not directly 
observable in the Im s correlation function; rather, its correlation length is 
given by the energy E which depends on the spatial momentum according to 
some dispersion relation. Exploiting the cubic symmetry of the Hamiltonian, 
which includes space inversion and rotations by it/2 about the coordinate 
axes, one finds that an inverse propagator for the Potts model must take the 
form m 2 + p\ + P2 +P3 + C(j° 4 )- Hence, identifying p% = iE, the exponential 
decay of the correlation function is related to the particle mass and the spatial 
momentum by 

E 2 = m 2 +pl+p 2 2 + 0(p\p 2 E 2 ,E 4 ) (2) 

In this Letter, eq. (|j) will be used only at minimal momentum: p$ min = or 
7r/24, depending on whether a periodic or C-periodic boundary condition is 
imposed in the corresponding direction. Furthermore, the energy levels we 
consider will be of order 0.1 in lattice units. Thus 0(p 4 ,p 2 E 2 , E A ) amounts 
to a correction in the 1% range. This is unresolvable within our statistical 
errors and will be neglected. 

Since the C-periodic boundary conditions allow only one ordered state 
for (3 > /3 C , vacuum tunnelling can take place only in the coexistence region 
of width o L around (3 C ^ on a finite L 3 lattice. As a consequence, possible 
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contamination of physical mass by the tunnelling mass is restricted to this 
region only. Figure [l] shows schematically the expected behaviour of the 
lowest excitation energy, near the phase transition, as a function of (3. Also 
shown is the corresponding behaviour for the periodic boundary case ||. 

Following Wang and DeTar we estimate the tunnelling effects for 
our version of the Potts model in terms of a simple model. Compared to 
their case, the Potts model with C-periodic boundary conditions has a much 
reduced symmetry — only charge conjugation symmetry is left. We, therefore, 
choose an extensive "ordered" state |0) which derives from the exact ordered 
vacuum at infinite (3, and another extensive state |D) which derives from the 
exact disordered vacuum at (3 = to descibe our model. Both these states 
are assumed to be invariant under charge conjugation. In case of |0) this 
is because all spins s are equal to the real number 1; in case of |D) it is 
due to the superposition of all spin eigenstates with equal amplitudes. We 
furthermore assume the following relations to hold for the matrix elements 
of the real part f of the spin operator: 

(0|r|0) = l (0|f|D) = (D|f|D) = 

When restricted to the Hilbert space spanned by |0) and |D) the Hamiltonian 
will take the form 

( E c \ 
[ c E D ) 

where c > as a matter of phase convention for the two states. At the point of 
order-disorder coexistence the energies Eo and Ed are equal. More generally, 
they depend on (3 and are proportional to the spatial volume V except for 
boundary effects. The matrix element c describes the vacuum tunnelling to 
be expected for V < oo; it will depend on (3 and will decrease with V — > oo. 
For the eigenstates of the Hamiltonian are superpositions of the form 

|+) = sina|0) + cosa|D) 
|— ) = cos a\0) — sin a|D) 

Eq — Et 



where tana = yl + x 2 + x with x = (3) 

2c 

Hence we have 

(+|r|+) = cos 2 a (— \r\— ) = sin 2 a (+|f|— ) = ( — |f|+) = sina cosa 
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The energy eigenvalues are 

E± = Eo + Ed ± £ where e = ^(E Q - £ D ) 2 + 4c 2 

Let |— ) correspond to the lower energy level which we put to zero, and let 
|+) be the first excited state at energy e. The partition function is 



1 + e 



-/3e 



Evaluating traces in the |±) basis we find for the statistical expectation value 
of the real part of the spin 



(r) 



cos 2 a + e /3e sin 2 a 
1 + 



To recover the canonical picture, let us assume that the spatial volume is 
finite but large. At the coexistence point, e ~ so that (f) « |. Elsewhere, 
e oc V, hence the exponential factors are very small, and (f) « cos 2 a. 
Furthermore, eq. (||) implies a ~ for E Q < E D (broken phase) and a ~ ~ for 
-Eo > Ej) (symmetric phase). Thus the spin expectation value indeed meets 
our expectations in this vacuum model. As for the correlation function, 

(f(O)r(r)) = Z- 1 Tr (fe~ TH re rH e'^ 11 ) 

we have 

N v, sin 2 a cos 2 a , . , „. cos 4 a + sin 4 ae _/3e 



Right at the phase transition where cos a = sin a = y 1/2 this expression 
simplifies to 

w«w \ \ 1 coshe(r - i/5) „ 

r(0 r(r = 777^7^ + t- 2 

4 cosh(/5e/2) 

As for particle excitations, we assume that one-particle states are generated 
by acting with the local spin operator s(x) on a vacuum state. The real 
part f(x) of the spin field operator has C-parity +1 while the imaginary part 
i(x) has C-parity — 1. The Hamiltonian H and the states |0) and |D) are 
charge-conjugation invariant. In the spin-spin correlations we thus expect to 
see an excitation scheme whose levels can be classified by their C-parity. 
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These model considerations suggest that vacuum tunnelling in the C- 
periodic Potts model (order-disorder tunnelling) only occurs in a narrow 
region around the phase transition, as illustrated in Figure [I]. Outside that 
region, both in the ordered and disordered phase, the energy spectrum con- 
sists of a nondegenerate extensive ground state, followed by particle excita- 
tions. In the simplified model described above, e is the energy of the first 
excited state which corresponds to the physical mass sufficiently away from 
(3 C . In the coexistence region, however, e ~ c and is thus governed entirely 
by the finite volume effects. Hence on increasingly larger volumes e becomes 
vanishingly small at the transition point. On the other hand, the width of 
the coexistence region also goes down linearly with volume, thus opening the 
possibility of obtaining the physical mass at (3 C in a limiting procedure. 

We have simulated the 3-dimensional 3-state Potts model on a 24 x 24 x 48 
lattice, imposing C-periodic boundary conditions in one of the spatial (L = 
24) directions, at couplings (3 = 0.5505, (3 = 0.551 and (3 = 0.55125. A total 
of 40000 measurements were performed at each coupling. At /3 — 0.5505 
and j3 = 0.551 measurements were interspersed with 20 updatings of the 
spin configuration; at /3 — 0.55125 with 100 updatings. Plots of the MC 
results for those couplings are very similar. We here present plots only for 
/3 = 0.5505. 

Figure |2| shows the Monte Carlo history of the real part of the spin av- 
eraged over the 3-dimensional lattice. The Figure illustrates the tunnelling 
between a globally ordered and a globally disordered state. 

Our MC data for the Res two-point correlation function R(d), pro- 
jected onto zero spatial momentum (Figure |3|) are strongly dominated by 
a constant — presumably the one modelled by eq. (|4]). The measurement- 
to- measurement fluctuations in R(d), also, are similar for all values of d, 
supporting the interpretation as a fluctuating constant. In fact, the errors 
in the exponentially decaying particle contributions which we are interested 
in are much smaller. To extract the latter, a better observable, suggested 
by eq. (|j), is the differences R(d) — R(d + 1). The constant part of the spin 
correlation thus drops out from both the mean values and the error bars, as 
is evident from Figure |j. A one-particle contribution to this difference would 
be of the form Asinhm(^Np — ~ — d) where Np is the temporal lattice size. 
Figure [| shows the corresponding fit of the data. The mass of the excitation 
is given in the Table. 

Practically the same mass is obtained by a global \ 2 fit °f the form 



8 



A cosh. m(Np — d) + B to the unsubtracted function R(d). In this latter 
way one obtains an excessively shallow x 2 minimum since x 2 is then mainly 
based on the fluctuations of the constant. Nevertheless, the agreement of 
mean excitation energies obtained from our R(d) and from R(d) — R(d + 1) 
supports the conjecture that C-periodic boundary conditions eliminate order- 
order tunnelling so that vacuum tunnelling effects must be taken into account 
only near the phase transition (order-disorder tunnelling). 

For the correlation of the imaginary part of the spin, projected onto 
minimal spatial momentum, Figure || shows the differences 1(d) — I(d + 1) 
and a fit of the form A sinh E(\Np — | — d). The mass of the C= —1 particle 
as obtained from formula (H) is given in the Table. In view of the non-zero 
momenta involved, we expect to find, as a consistency check, no particle 
energies in the Im s correlation 

(a) below tt/L if we impose C-periodicity in only one spatial direction, 

(b) below \j2it 2 / L 2 if both spatial directions are C-periodic. 

In our sample computations of case (a), L = 24 so that the physicality 
bound is 0.131. This bound is consistent with the x 2 distributions we obtain 
at all couplings from the differences 1(d) — I(d +1). In the unsubtracted 
correlation functions 1(d) we actually observe small constant contributions 
as well. Furthermore, at all three couplings we find error bars consistent with 
the fluctuations of the mean values only in 1(d) — I(d + 1) while for 1(d) the 
mean values are overlaid again with the error bars of a fluctuating constant. 
However, the size of the effect is such that it can easily be understood in terms 
of mere statistical fluctuations. The constant fitted to 1(d) at (3 = 0.5505, 
for example, is 1.2 x 10 -4 which is less than a standard deviation, after 
40000 measurements, of a fluctuating constant with zero average and an 
RMS of the same order as in R(d). In the following, we shall form differences 
1(d) — I(d + 1) with the intention to eliminate from 1(d) the artifacts of a 
finite number of iterations. 

At (5 — 0.551 and (3 = 0.55125, where we expect to be sufficiently deep 
in the ordered phase, the maximum-likelihood estimates (mean values in 
Table 1) of m_ = 0.11 for both cases indicate a massive particle in the 
C = — 1 channel degenerate with the C = +1 particle discussed above. 
Since according to formula (0) involves a difference, E 2 —p 2 min , the mass 
estimates in the C = — 1 channel can be quite sensitive to uncertainties in the 
correlation length £ = E^ 1 . At the above couplings this results in a rather 
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large error estimate for m_. At /3 = 0.5505, i.e. near the phase transition, 
we find an essentially unchanged m_ = 0.13 ± 0.02, while m + decreases to 
0.07 indicating the increasing dominance of order- disorder tunnelling in this 
regime. 

In conclusion, we have shown that the C-periodic boundary conditions, 
given by eq. ([!]), eliminate order-order tunnellings and allow numerical im- 
plementation of the standard statistical-mechanical definition of connected 
correlation functions and particle masses in the 3-dimensional 3-state Potts 
model. This is also true for all Z(N) models, A-state Potts models and 
SU(N) gauge theories at finite temperature in any dimensions, provided iV 
is odd. In particular, quenched QCD also falls in this category. Further- 
more, we suggest a method to obtain the physical mass in all these theories 
at (3 C in the thermodynamic limit. It will be interesting to find out using this 
method whether the first order phase transition in the 3-dimensional 3-state 
Potts model or the 3-dimensional SU(3) gauge theory at nonzero tempera- 
ture manifests itself as a discontinuity in the physical correlation length. 

Also, it will be interesting to analyse the dispersion relation (|2|) beyond 
the quadratic order in the energy and momenta. For quasi-free actions with 
a Laplacian structure of derivatives, E 2 = m 2 + J2k sm2 Vk is a natural 
extension; it has, in fact, produced satisfactory mass corrections for U(l) 
monopoles in the Villain approximation on smaller lattices ||. For the Potts 
model, the action is of an entirely different form. Here one can study nu- 
merically the correlation functions of the real and imaginary parts of the 
spin after projection onto various non-minimal spatial momenta. However, 
on lattices large enough so that 0(p^ in ) i s negligible, the Ims correlation 
function in conjunction with dispersion relation (|2]) offers perhaps the best 
possible way to obtain the mass of the lowest C — — 1 particle as it is not 
overlaid with any vacuum tunnelling. 
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Figure 1: Schematic diagram illustrating the behaviour of the energy of the 
first excited state as a function of f3 near the phase transition for two volumes 
Vi < V 2 , using C-periodic (full lines) and periodic (full and dashed lines) 
boundary conditions. 



Figure 2: Monte Carlo history of the real part of the spin (3-dimensional 
average) at (3 = 0.5505. 



Figure 3: Monte Carlo data for the Res correlation R(d) at (3 = 0.5505 and 
fit of the form Acoshm(24 — d) + B. The error bars are dominated by the 
fluctuations of the constant B. 



Figure 4: Monte Carlo data for the differences R(d) — R(d+ 1) at f3 — 0.5505 
and fit of the form Asinhm(23.5 — d). The x 2 fit is based on data points 
5 < d < 23. 



Figure 5: Monte Carlo data for the differences 1(d) — I(d+ 1) at /3 — 0.5505 
and fit of the form Asinhi?(23.5 — d). The x 2 At is based on data points 
5 < d < 23. 
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Table 1: Inverse correlation lenghts £~ 1 and particle masses obtained from the 
spin-spin correlation. R and I refer to the correlation of the real or imaginary 
parts of the spin, respectively. The form of the fit is Asinh^ _1 (23.5 — d). 



p 


Function 


r 1 


Particle Mass 


0.55050 


R(d) - R(d + 1) 


072 +0 - 012 
UAJ ^ -0.014 


same as 


0.55050 


I(d)-I(d+1) 


nio , +0.016 
°' 184 -0.014 


0.13 + 0.02 


0.55100 


R(d) - R(d + 1) 


108 +0 -° 12 
°- 108 -0.014 


same as 


0.55100 


I(d)-I(d+1) 


17 +0 -° 5 
-0.045 


on + - 07 

-0.11 


0.55125 


R(d) - R(d + 1) 


134 +0 -° 08 
-0.007 


same as 


0.55125 


I(d)-I(d+1) 


01? +0.045 
-0.055 


+ 0.06 
-0.11 
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